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ABSTRACT 

The rapidly declining population of bright quasars at z J 3 appears to make an increasingly small contribution to the ionising 
background at the H i Lyman limit. It is then generally though that massive stars in (pre-)galactic systems may provide the additional 
ionising flux needed to complete H i reionisation by z > 6. A galaxy dominated background, however, may require that the escape 
fraction of Lyman continuum radiation from high redshift galaxies is as high as 10%, a value somewhat at odds with (admittedly 
scarce) observational constraints. High escape fractions from dwarf galaxies have been advocated, or, alternatively, a so-far undetected 
(or barely detected) population of unobscured, high-redshift faint AGNs. Here we question the latter hypothesis, and show that such 
sources, to be consistent with the measured level of the unresolved X-ray background at z = 0, can provide a fraction of the H n filling 
factor not larger than 13% by z — 6. The fraction rises to < 27% in the somewhat extreme case of a constant comoving redshift 
evolution of the AGN emissivity. This still calls for a mean escape fraction of ionising photons from high-; galaxies > 10%. 
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1. Introduction 

The reionisation of the all-pervading intergalactic medium 
(IGM) is a landmark event in the history of the Universe. Studies 
of the so-called Gunn-Peterson absorption in the spectra of dis¬ 
tant quasars show that h ydrogen was already highly ioni sed out 
to redshift z ~ 6 (e.g., Songaila '2004b iFan et al.l120061) . while 
CMB polarisation data constrain the redshift of a sudden reion- 
isation event to be signifi cantly higher, z ~ 10 dJarosik et al.l 
1201 ltlHinshaw et al.ll2013l) . 

Most of our understanding of IGM physics, and its im¬ 
plication for galaxy formation and metal enrichment, depends 
critically on the properties of the cosmic ionising background. 
While it is generally thought that the gas is kept ionised 
by the integrated UV emission fro m ac tive galactic nuclei 
(AGNs) and star-forming ga laxies (iMiralda-Escude & Ostriken 
ll990i:lHaardt & Madaull 996i) . the relative contributions of these 
sources as a function of cosmic time are poorly known. 

At z > 3, the declining population of bright quasars ap¬ 
pears to make an increasingly small contribution to the ion¬ 
ising radiation background at the H i Lyman limit. It was 
then suggested that massive stars in galactic systems may 
provi de the additional ionising flux needed at early times 


(e.g. Madau et al. 1999; Gnedin 

2000; Wvithe & Loeb 

2003; 

Meiksin 2005; Trac & Cen 2007 

Faucher-Giguere et al 

2008; 

Gilmore et al. 2009; Robertson et al. 20IQ). However, leaking 


Lyman continuum radiation from bright galaxies seem to be 
modest (see, e.g.. IVanzella et all 120101) . and it has been there¬ 
fore argued that dwarf galaxies (with virial mass below ~ 10 9 
M 0 ) may produce the dominant contribution t o the H i ionising 
UV background (e.g jRobertson & Ellisll2012b . 

Alternative to invoking a major contribution to reionisa¬ 
tion from dwarf(ish) galaxies is the possibility that the AGN 


emissivity at z > 4 is indeed much larger than generally 
thought. Ind ications along such line have been reported by sev¬ 
eral g roups dGlikman et all201lLlCivano et alJl201 lUFiore et al.l 
12012l) . thou gh it is fair to say t hat results see m not so univo¬ 
cal (se e, e.g.. (Masters et al.l2012l) . Very recently, Giallong o et al.l 
(12015h found 22 AGN candidates at z > 4 in the Candel/GOOD- 
S/Chandra Deep Field South field, suggestive of a prominent 
contribution of AGNs to the ionising background in the range 
4 < z < 6.5. The resulting H i photoionisation rate is in¬ 
deed consistent with various estimates at the same redshifts, 
based on both the flux-decrement and proxim ity effect tech¬ 
niques (iBecker et all2007t [Calverlev et al.ll201 lh . 

The high redshift population of AGNs should leave an im¬ 
print in the observed cosmic X-ray background (XRB). Chandra 
deep observations resolved the XRB into discrete sources at a 
level of 80 - 90% over the entire bandwidth (0.5 - 2 keV), 
with only a fra ction ~ 1% of the signal arising fro m sources 
located at z > 4 (IXue et al.ll201 lIlTMoretti et al.l(l2012l) exploited 
the very low instrumental noise of the Swift XRT to measure 
thejrtill unresolved XRB spectrum at the highest accuracy. In 
ISalvaterra et al.l (20121 ) we used such measures to place upper 
limits on the cosmic accretion history of massive black holes. 
An obvious caveat to our conclusions is the possible existence of 
a large population of severely obscured (log Nh f 25)accreting 
black holes, hence not glowing in the X-rays. Unless advocating 
a very peculiar UV-to-X-ray spectral energy distribution, such 
caveat would not apply if the unresolved AGN population does 
contribute significantly to the ionisation background. The as¬ 
sessment of the contribution to the XRB of such UV-emitting 
AGNs is precisely the goal of this Letter. Specifically, we will 
translate the lMoretti et al.l (12012i) upper limits to the unresolved 
XRB into upper limits on the possible contribution of high- 
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redshift, unobscured fain t AGNs to H i reioni s ation. A similar 
analys i s was propo sed by Diikstra et al.1 (120041) . ISalvaterra et alJ 
(I2005L 120071) . and McOuinnl (1201 2l) . with conflicting results. 
Here we use the most updated limits on the XRB and adopt a 
(h, Cl m , Qa) = (0.7,0.3,0.7) cosmology. 


2. Methodology 

Assuming an AGN comoving X-ray specific emissivity oc 
(1 + zY y , the XRB at observed energy Eq due to sources 
located at redshift z> z x is: 


Je 0 (> Zx) 


J r»cx 


dz 


ftkeV 


f—V 

\ El keV / 


(l +zy 


( 1 ) 


where E = E f) ( 1 + z). We now relate the specific emissivity at 
2 keV, t 2 kev 5 to that at 912A, i,.e., e 2 keV = Ke 912 ^, where the 
“/("-correction" normalisation reads 

_ / 1300A f uv /2500A\ ffuv / £ 2 k e v \~ am ... 

\ 912A / \ 1300A/ U2500A/ ' 

Here a ox is the optical-to-X-rays spectral index, defined from 
the specific emissivity at 2 keV and at 2500 A, a ox = 
-0.384 log(e 2 kev/f 250 oA)- I n writing eq.[2]we followed the piece- 
wise UV AGN spectr al energy distribution as described in 
lHaardt & Madau l ( 2012 ). Now the r.h.s. integral in eq.|T|is easily 
solved, and the obtained Je 0 (> z x ) constrained to be not larger 
than the observational upper limit /j* s . This in turn gives the 
maximum value of e gl2 A consistent with the limits on the unre¬ 
solved XRB: 

^ 77(1+ + zyy (3) 


where we neglected the energy density of the cosmological con¬ 
stant (we are interested in the redshift regime z » z m A — 0.33). 
We set 1 ; = (y + a x + 3/2) and the term Ru > 1 is meant to 
account for the contribution of obscured AGNs at z > Zx to the 
XRB observed at energy E g . 

It is now straightforward to translate the above limit into a 
limit on reionisation. The volume filling factor of H 11 regions 
is the solution of the following differential equation (see 
Mad auet alii 19991) : 


dQmi _ n_ _ Qvat 

dt tl 0 free 


(4) 


where no is the cosmic hydrogen mean density, and n(z) — 
e 9 \ 2 A(z)l(h p af m ) the photon emission rate ( h p is the Planck con¬ 
stant, and the FUV emissivity is oc y~ Q,fllt ). The H 11 re combina- 
tion time f rec is computed as in lHaardt & Madaul(i2012l) . 


3. Observational Parameters 


The upper limits given by eq.[3]and eq.[4]depend upon a number 
of parameters which need to be observationally constrained. In 
this section we discuss the choice we make for each of them. 

Unresolved XRB. The unresolved XRB shows a very hard 
spectrum, suggesting that most (if not all) of the flux comes from 
low-z obscured AGNs. Indeed, Moretti e t aL ( 20121) . by adopting 
the XRB synthesis model of iGilli et all ( 2007 ). derived a strin¬ 
gent limit /° bs < 1.9 x 10- 27 erg cm -2 s~ ] Hz" 1 sr" 1 at Eq - 1.5 
keV once accounting for absorbed AGNs at z IS 5 whose fluxes 
lie below the Chandra limit. Yet, the synthesis model falls short 


at Eo > 3 keV, suggesting the possible existence of a population 
of Compton thick AGNs at intermediate redshifts. 


UV and X-ray spectral ind ices. As already stated, we use the 
very same parametrisation of lHaardt & Madaul (120121) . adopting 
for the UV slope a^v = 0.44 (T > 1300 A. IVanden Berk et al.1 
(2001 )), an d af uv = 1.57 in the FUV range (A < 1300 A, 
iTelfer et ak d2002l) ). Concerning the X-ray spectrum, we no¬ 
tice that most of the contribution to the XRB at £0 = 1.5 keV 
is expected from source located just above z — 5, i.e., from 
photons emitted at rest-frame energy - 10 keV. In this energy 
range unobscured AGNs exhibit a power-law spectrum with in¬ 
dex a x - 0.8, as a combination of the intrin sic c ontinuum and 
of the Compton reflection bump (see, e.g.. llJeda et all ( 20141) ). 
In our analysis we then take a x = 0.8. From eq. [2] it is apparent 
how the exact values of the UV and X-ray spectral indices will 
affect our conclusions only marginally. 


Optical-to-X-ray spectral index. The value of a ox has a ma¬ 
jor impact on our estimate of Qhii since (LAkev/lQsooA) — 403. 
The study of the correlation between X-ray and UV luminosities 
has been the subject of many works on both optically selected 
and X-ray selected AGNs. Among others, ISteffen et alJ (120061) 
found a significant correlation between a m and the monochro¬ 
matic luminosities at 2500 A (i.e., Ti 500 a L 2keV , with ft > 1) 
in a sample of 333 optically selected AGNs. No significant cor¬ 
relation of a ox with redshift was reported. iLusso et aD ( 1201 O t) 
analysed a sample of 545 X-ray selected Type I AGNs from 
the XMM-COSMOS survey, finding again a correlation between 
a ox and L 2S00A . The mean value of a ox for the full samp le was 
1.37 with a dispersion around the mean of 0.18. iMarchese et aD 
( 2012 1) also reported a highly significant correlation between a ox 
and the UV luminosity in a sample of 195 X-ray selected Type 
I bright AGNs, basically confirming the ILusso et alJ ( 2010 ) re¬ 
sults. In our investigation the supposedly unaccounted popu¬ 
lation of AGNs responsible of H 1 reionisation must necessar¬ 
ily resides in the very faint-end of the UV luminosity function. 
According to the literature cited above, this would imply a value 
of a ox on the lowest side of the distribution, though it must be 
considered that the redshift range of interest here is basically un¬ 
explored at X-ray wavelengths. Given that, we adopt a fiducial 
value a ox = 1.35. We are confident that, if the observed corre¬ 
lation between T 2 50 oA anc * a ox holds at very high redshifts, such 
choice is conservative, hence strengthening our conclusions. 


Redshift evolution. The evolution of the AGN space density 
at high redshifts has been subject of several revisions in the last 
decade, mainly because of the dearth of data at z > 4. As an 
example. llJeda et akl ( 2003 ) adopted, for very luminous sources, 
an evolution factor oc (1 + z)~ y with y = 1.5 above z — 1.9, 
while fo r fainter AGNs the turn o ver occurs at increasingly lower 
redshift. Silverman et al . (200 8) found a much sharper decline, 
y = 3.27, similar to that derived in studies of optically selected 
QSOs. Recently. iHiroi et aTl (t20 12 ) claimed an even stronger de- 
cline at z > 3, y — 6.2, a value adopted bv llJeda et alJ (120141) in 
the most recent and updated study of the hard X-ray LF. The 
situation is somewhat more confusing in the optical-UV band. 
While different groups agreed on the faint-end slope of the LF, 
- 1.7, they sorted out quite different absolute space densities. 
Specifically, iGlikman et al.1 (1201 ll) cl aim roughly a fa ctor four 
more so urces at z > 3 comp ared to [ikeda et alJ (1201 ll) . More 
recently. Masters et Id] (1 201 2l) found a decrease by a factor of 
four in the number density of faint QSOs in COSMOS between 
z ~ 3.2 and z ~ 4, suppo rting the result s of Ikeda et aTl (1201 ll) . 
Overall, the results from [Masters et al.l (1201 2l) suggest a simi¬ 
lar evolution of the UV and X-ray LFs at z > 3. However, 
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a large normalisation of the UV LF, basically consistent with 
iGlikman et al.1 (1201 ll) though at higher redshifts (z — 5 - 6) and 
fainte r UV magnitudes, was recently claimed (Giallong o et alJ 
120151 ). Given these uncertainties, we will assume a redshift evo¬ 
lution of the emissivity as sharp as y — 6.2 in both the X-ray 
and UV bands, but we will also show results for the somewhat 
extreme case of a constant comoving emissivity (y = 0). 

Obscured sources. The parameter Ru in eq. [3] is meant to 
account for the contribution to the XRB given by sources ob¬ 
scured in the optical band, thus not contributing to the ionising 
background. Such contribution could be relevant since photons 
observed at Eq = 1.5 keV are emitted, by z > 5 AGNs, at rest 
frame energies > 10 keV, where the emission of absorbed AGNs 
is anyway relevant. We implemented the X-ray LF of Ueda et al J 
(120141) . and found that JeJ> z x ) (with Eq = 1.5 keV and z x = 5) 
is almost evenly divided between objects with log Nh < 22 and 
objects with log Nh > 22, which would give Ru - 2. However, 
it is not simple to determine above which X-ray determined 
equivalent hydrogen column density Nh sources are severely 
obscured in the optical-UV band. In their study, [Masters et all 
(120121) found that ~ 75% of X-ray bright AGNs at z ~ 3 - 4 are 
indeed optically obscured. Taken at face value, this would imply 
Ru — 4. Finally it is worth noticing that at lower redshifts (z < 3) 
the incidence of obscured AGNs is strongly anti-correlated with 
X-ray luminosity dMerloni et al.ll2014 ). Provided that the trend 
is similar at earlier epochs, this very fact points toward a high Ru 
correction factor. In our fiducial model we then assume Ru = 4. 

Lower redshift of unresolved XRB. In our analysis, the limit¬ 
ing redshift z x plays an important role, as it sets the minimum 
redshift of the unaccounted AGN population we are testing. 
Clearly, such population must give a contribution to the XRB 
not exceeding the measured unresolved fraction. Specifically, 
the upper limits to the unresolved XRB given by IMoretti et all 
d2012l) were obtained subtracting to the total XRB all sources 
listed in the 4Ms-Chandra catalog (IXue et al.ll201 ll) . which basi¬ 
cally contains no AGNs with z > 5. As and e xamp le, among the 
6 AGN candidate s at z > 5 found by iGiallongo et all (120151) . 
only 2 are in the IXue et all (1201 ll) catalog. Should the unre¬ 
solved XRB arise from high-z AGNs, it must necessary come 
from z > 5 sources unless they own a very peculiar redshift dis¬ 
tribution. Given that, we assume z x = 5 as the lower limiting 
redshift in our study. 

To summarise, our benchmark model adopts: (T| UV = 1.57, 
a uv = 0.44, a x = 0.8, a ox = 1.35, y = 6.2, R n =4, z x = 5, and 
the XRB limit 7°j S keV < 1.9 x 10 27 erg cm 2 s 1 Hz 1 sr 1 . 


4. Results 


Fig. Q] shows e 912 A given by eq. [3] The emissivity at the Lyman 
limit is a most interesting quantity, as it can be compared to di¬ 
verse observational estimates existing in literature. Our bench¬ 
mark case is shown as the red solid line starting from z = 5. Our 
limit is compared to the recent values of Giallongo et all (120151) 
and to the values reported by [Masters et al.i d2012l) (shown as 
black data points and open triangles, respectively). An assess¬ 
ment of the reasons behind the discrepancy between these dif¬ 
ferent results is beyond the scope of this letter. Still, taking the 
[Giallongo et al.l t2 015l ) ionising emissivity at face value, we must 
conclude that the associated AGNs basically saturate the ob¬ 
served XRB, as apparent from Fig. |T] 

We also compare our estimate of f,, , A with the “minimum 


reionisation model” of Haardt & Madau (120121). In Fig. Q]the 
overall Lyman limit emissivity of iHaardt & Madau (120121) is 



Fig. 1 . The maximum emissivity at the H 1 Lyman limit vs. red¬ 
shift. The red solid curve at z > 5 is our benchmark case, consis¬ 
tent with the XRB limit 7°5 S keV < 1.9 x 1 0 -27 erg cnT 2 s 1 Hz -1 
sr -1 . Such value accounts for source s at z < 5 below the Chandra 
flux limit, estimated from the iGilli et akl (2 0071 ) XRB synthesis 
model. The dashed line shows the extreme case of a constant 
comoving redshift evolution of the AGN emissivity. The black 
line is the e missi vity employed in the UV background model of 
iHaardt & Madau ( 201 2). with galaxies (dotted green line) and 
AGNs (solid green line) shown separately. The black dat a point s 
are the AGN emissivity as estimated by GiallongoetalJ 00151) . 
while open triangles are data from lMasters et ali i 20121) . 


shown as a black solid line, along with the separate contri¬ 
bution of AGNs (green solid line) and star forming galaxies 
(gr een dotted line ) . The AGN emissivity clos ely fits the results 
by Hopkins et all (120071) . while for galaxies |Haardt & Madaul 
(2012) assumed that the fraction of Lyman continuum photons 
leaking into the IGM is a strong increasing function of redshift. 
Such model reionises Hi by z =* 6.7 and He 11 by z =* 2.8. 
Though o ur estimate benchmark e missivity at z - 5 - 6 is sim¬ 
ilar to the [Haardt & Madaul d2012l) model, the approximatively 
constant comoving behaviour of the latter compared to the steep 
decline we adopt here leads to the different outcome in terms 
of reionisation (see next). For such reason we tested the some¬ 
what extreme case of a constant comoving emissivity, shown as 
a dashed line in Fig. D] The re s ulting maximum emissivity is 
lower than the [Haardt & Madaul (120121) reionisation model by a 
factor of =s 4 at high-z. 

In Fig. [2] we show the resulting volume fraction occupied 
by H 11 regions (eq. [4]). Our benchmark case allows for a fraction 
< 13% of the IGM to be ionised by AGNs by z > 6, showing that 
AGNs alone can not reionise the Universe. We checked that this 
conclusion holds in spite of the uncertainties of the parameters. 
Qhu = 1 at z = 6 can be reached only for a ox > 1.7, or a ox > 
1.5 assuming Ru — 1 (i.e., no obscured sources). Such figures 
seem to be unlikely when compared to available data. Finally, the 
constant comoving case (dashed line) produces a more extended 
reionisation history, still it can only account for < 27% of the 
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Fig. 2. Maximum volume filling factor of H u vs. redshift. 
Curves as in Fig. Q] The dashed vertical line marks a fiducial 
reionisation redshift, z = 6. In our benchmark case the contribu¬ 
tion to H ii reionisation from AGNs is below 13%. 


ionised volume at z > 6. In this case, reionisation by z > 6 
requires a mean escape fraction from star forming galaxies > 
10 %. 

It is interesting to note that_the benchmark case (as well 
as the AGNs observed by iGiallongo et alJ ( 2015 1) would pro- 
duce a H i ionisation rate consistent with z — 5-6 data (e.g., 
IWvithe & Boltonll201 it ICalverlev et alJl20llh . still it falls short 
in reionsing the IGM at z > 6 (Fig.0. In other words, matching 
the observed level of the ionising background just below the ion¬ 
isation redshift does not guarantee that a particular model is ac¬ 
tually consistent with the entire reionisation history of the IGM. 

A final comments concerns He n reionisation. Though a de¬ 
tailed assessment of such process is beyond the scope of this 
Letter, an AGN dominated background would certainly lead to 
an extended reionisation epoch. This may agree with the re¬ 
cent claiming of IWorseck et ali (i2014l) . but it may be in conflict 
with the sharp increase of the IGM temperature at mean cosmic 


density observed in the range 2 < z < 4 (Schave et a 

j 2000; 

Becker et al. 2011: Bolton et al. 20121 2014; Boera et al. 

2014. 


5. Conclusions 

Under reasonable assumptions, we have shown that a population 
of unobscured, UV emitting AGNs at z > 5 if leading H i reion¬ 
isation would exceed observational constraints derived from the 
unresolved fraction of the X-ray background. Even a constant 
comoving emissivity at high-z would not be enough to produce 
an AGN dominated ionising background. AGNs can account for 
a fraction of the ionising photon budget < 13%, calling for a 
dominant contribution from star forming galaxies. Given the ob¬ 
servational constraints (e.g jBouwens et akil201 11 1 on the galaxy 
population at high-z, this in turns requires a large mean escape 
fraction (> 10%). 
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